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We consider small ballistic quantum dots weakly coupled to the leads in the chaotic regime and 
look for significant spin-orbit effects. We find that these effects can become quite prominent in the 
vicinity of degeneracies of many-body energies. We illustrate the idea by considering a case where 
the intrinsic exchange term — JS 2 brings singlet and triplet many-body states near each other, 
while an externally tunable Zeeman term then closes the gap between the singlet and the one of the 
triplet states (with spin projection parallel the external field). Near this degeneracy, the spin-orbit 
coupling leads to a striking temperature dependence of the conductance, with observable effects of 
order unity at temperatures lower than the strength of the spin-orbit coupling. Under favorable 
circumstances, spelled out in the paper, these order unity effects in the conductance persist to 
temperatures much higher than the spin-orbit coupling strength. Our conclusions are unaffected 
by the presence of non-universal perturbations. We suggest a class of experiments to explore this 
regime. 
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I. INTRODUCTION 



Electrons in planar quantum dots experience spin- 
orbit interactions due to intrinsic (Dresselhaus 1 ) terms 
and due to the electric field arising from the confinement 
to two dimensions (the Rashba term 2 ). Their strength is 
measured by the parameter L/X, where L is the dot size 
and A is a spin-orbit scattering length. (We do not worry 
in this paper about the two different scattering lengths 
for the two terms.) Spin-orbit couplings have been thor- 
oughly investigated in noninteracting quantum dots 3-9 , 
and more recently, in interacting dots as well in certain 
special regimes 10-12 . Spin-orbit terms violate spin con- 
servation and appear to move the symmetry class from 
orthogonal to symplectic to first order in L/X. Aleiner 
and Fal'ko (AF) 4 used a judicious unitary transforma- 
tion to show that despite appearances the spin-orbit in- 
teraction enters the hamiltonian only to order (L/X) 2 
and beyond. They also showed that to order (L/X) 2 (but 
not beyond) a new quantity af F (similar algebraically 
to the usual o~ z but distinct from it) is conserved. The 
crossover to the new symmetry class occurs at a scale 
e AF ~ (L/X) 4 gS, and the violation oiaf F occurs only at 
the scale 4 (L/X) 6 g5, where 6 is the mean single-particle 
level spacing, g5 is the Thouless energy, and g is the 
Thouless number. Thus in a small dot with L ~ 200 nm 
and A ~ 2000 — 3000 nm, spin-orbit effects are expected 
to be extremely small. 

Here we ask how these considerations, derived for non- 
interacting particles, are modified if the exchange inter- 
action (— JS 2 ) of the Universal hamiltonian 14,15 is in- 
cluded. Recall that the Universal Hamiltonian is the 



correct low-energy fixed point 16 deep inside the Thouless 
band for sufficiently weak interactions 17 . We find that in 
the presence of the exchange coupling J the spin-orbit 
effects really are of order in the hamiltonian and 

that in the proximity of many-body degeneracies (in- 
duced largely by J and aided by a Zeeman field which 
also serves as a knob for dialing through the transition), 
effects linear in L/X may be found in energies, and most 
strikingly, effects of order unity may generically be found 
in the conductances (provided T is low enough). This 
paper focuses on calculating these effects, and proposing 
experiments to measure them. The optimal situation is 
in dots where low temperatures (in units of level spac- 
ing) may be attained without too much reduction of 4. 
Currently, for dots of linear size L ~ 200nm, the typical 
mean level spacing is 5 ~ 70— lOOueV and temperatures 
of T ~ .1(5 are achievable. The effects under study will 
be more readily detected if T/ 5 is lower by a factor of 5 



In Section II we consider the Universal Hamiltonian 
plus a spin-orbit term. We review some previous work 
and elucidate how and when various powers of L/X ap- 
pear in observables and how the presence of a J ^ 
affects this. We then turn to our main result, which 
is to show that in certain conductance measurements 
there is a good chance of seeing spin-orbit effects of or- 
der unity. We estimate the effects in terms of a single 
parameter and propose experiments to reveal them and 
to interpret the data. Discussions and conclusions follow 
in Section III. 



1 



II LOW TEMPERATURE PHYSICS 



2 



II. LOW TEMPERATURE PHYSICS 



Our starting point is the hamiltonian 



H — Hjj + H so + H 



(1) 



The Universal Hamiltonian 14 for the case of orthogo- 
nal symmetry is 



(2) 



Here s creates an electron in a single-particle state \x 
of energy and spin projection s z = s/2 in a chaotic but 
ballistic dot. The charging energy is U, N is the particle 
number, and S the total spin. The essence of the spin- 
orbit hamiltonian is captured in first-quantization by the 
expression 



H.„ = 



z • er x p 

2mA 



(3) 



where z is a unit vector along the z-axis. For simplicity of 
discussion we keep the Rashba 2 term only. The Zeeman 
term due to a parallel field (chosen to be along the x- 
axis) is Hz = —EzS x . 

We will be dealing with the case when the spin-orbit 
coupling is very weak, in the sense that the Aleiner- 
Fal'ko crossover scale (j) 4 gS <C 5, which is applicable 
to small dots. The opposite regime for large dots leads 
to quite different physics 11,12 , which we will comment on 
in the conclusion. 

Before plunging into details we outline our strategy. 
Imagine that there is no spin-orbit coupling, that we 
are at very low temperatures T ~ 0, and are sitting 
at the Coulomb Blockade 18 conductance maximum ob- 
tained by tuning the gate voltage V g to make an TV-body 
ground state with (S = 0, S x = 0) = (0, 0) degenerate 
with an N + 1- body ground state with (S = ^,S X = 
i) = {\,\)- Now imagine turning on the Zeeman term 
Ez to cause a ground state transition in the A-body 
sector so that (S = 1,S X = 1) = (1,1) beats (0,0). 
The — JS 2 term lowers E z , the point where all three 
states-(0, 0), (1, 1), {\, \)- become degenerate, by bring- 
ing down all three triplet states. At E* z the conduc- 
tance will have a different value as compared to cither 
side where only two states are degenerate. For example 
(as will be seen in detail below) in the case where all 
tunnelling widths are equal at both point contacts, the 
conductance at E* z will be | of its value of on cither side. 

How does a spin-orbit interaction modify this expec- 
tation? By eliminating the level crossing, and yielding 
a unique ground state in the iV-body sector, it will, in 
the illustrative example being considered, kill the bump 
entirely. This is an example of the order unity effect 
on the conductance due the spin-orbit interaction. If 
we now heat up the system so T ~ a where a de- 
fines the scale of the spin-orbit interaction, the bump in 



conductance, rounded by the elevated temperature, will 
reemerge, since the two states that avoided each other 
will contribute more or less equally. When we get down 
to details we will find that while the above mentioned 
behavior can indeed occur (provided all tunnelling rates 
are equal at both leads), so can many others, some even 
more dramatic, the actual behavior being sensitively de- 
pendent on the wave functions at the leads. How then are 
we to predict what to expect? We shall see that with a 
few additional conductance measurements one can ob- 
tain the relevant information on these wave functions 
and know in advance what to expect of a given conduc- 
tance peak in a given dot. 



A. Essential background 

Here we review some recent results on the powers of 
the spin-orbit coupling that should appear in noninter- 
acting dots. 

Let us begin with a simple non-interacting hamilto- 
nian: 

where V is the confining potential. The spin-orbit term, 
considered as a perturbation, is of order 1/A and we 
consider its impact on single-particle levels described by 
time-reversal invariant real wave functions. We will refer 
to powers of j as powers of L/X since in the end the 
L in the numerator will invariably appear in physical 
quantities. For future reference, we define the orthogonal 
basis of H 



(5) 



Going over to second-quantized notation we can ex- 
press Hq in its own orthogonal eigenbasis in the many- 
body Fock space as 



(6) 



and the total spin operator as S = | ^ s a ss i c^i . 

Because H so has no diagonal matrix elements in the 
basis of the unperturbed hamiltonian, we expect this 
term to appear to order (y ) 2 in the energy shift of levels. 
It was pointed out by AF that this is illusory 4 . Upon 
completing squares they bring the kinetic term to a form 



(P + ^ff) 2 



1 



2m 



4mA 2 



(7) 



when it becomes apparent that the particle is coupled to 
a spatially constant nonabelian vector potential. Using 
the gauge transformation 
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U = exp 



iz ■ r x er 



2A 



(8) 



and expanding it in inverse powers of y, they transform 
to new variables in terms of which the kinetic term be- 



comes 



where 



Haf — 



{PAF 



aj 2 



2m 



r x z 



4mA 2 



2A 2 



(9) 



(10) 



Terms of order (y) 3 and higher are not shown. 

Note that a z F and paf are not the same as a z and p, 
while r remains untransformed because it commutes with 
U. We can rewrite Haf as an unperturbed orthogonal 
Hamiltonian Hq,af (the first two terms in the following 
equation) perturbed by spin-orbit terms 



Haf = ^ + V(v) + H soA f 
2m 



(11) 



Going over to second quantization, we express Hq^af in 
its orthogonal eigenbasis 



Haf = ^ e-n, 



AF<$n t tdn,t + H, 



so.AF 



(12) 



where we have denoted the new fermionic operators d, d) 
to emphasize their difference from the earlier c, cr's. We 
now define a total spin operator in the H ,af eigenbasis 

as Saf — \ du tO'tt'dftj'- Clearly, Saf is not the same 
operator as S. In fact, 



S = U^SafU = S 



AF 



SSai 



(13) 



which shows that Saf is some complicated combination 
of the original spin and orbital degrees of freedom. 

Note that the vector potential is not fully gauged 
away. This must be so since a spatially constant 
nonabclian vector potential can have a nonzero field 
strength, which in our case is f^ v = [A^ , A„] ~ a z /X 2 . 
This is why in the end the particle is coupled to a mag- 
netic field of size (y) 2 and a sign which depends on s z . 
This coupling will contribute to energies or rates to or- 
der ( j) 4 . It is also clear that a z in the new AF basis (to 
be referred to as <J Z F ) is conserved by the aj_ term and 
violated only by higher order terms. 

It is instructive to ask how this correct physics would 
manifest itself if one did not make the AF transformation 
and worked perturbatively in the original orthogonal ba- 
sis. The heart of the equivalence between the two bases 
is the operator identity 



^ = -AHoM 

m h 



(14) 



which means that 



(/*|P|"> = "^-(C/i 



tv)W\v) 



(15) 



This identity was emphasized and exploited by 
Halperin et al 3 in their analysis of the interplay between 
spin-orbit and Zeeman couplings. To see explicitly how 
this enters our discussion, consider energy level shifts of 
order (y) 2 in the orthogonal basis. One finds for the 
second order energy shift in the single-particle state \i of 
either spin: 



AE^ = 1 
" 4m 2 A 2 



E 



MpH • (HpIm) 



4m 2 A 2 h 
1 



E 



( M |pH • (H [H , r] \ V ) 



(16) 
(17) 



im 



4m 2 A 2 h 
1 im 

~4m 2 A 2 ~?T 



((M| P T| M } + ( M |r| M }.( M |p| M }}) (18) 



(mIp-tIm) 



(19) 



where we have used the time-reversal* invariance of the 
orthogonal states |/x) to set (/x|p|/x) = 0. Using time- 
reversal invariance again we get 



(//|p-r|/x) =-(n\r-p\n) 



1 



( M |[p,-r]|/i) (20) 



ending up with a shift — 4 J; A ^ which is the state- 
independent shift obtained by completing squares. 

At fourth-order (-^) 4 one will find corrections that are 
state -dependent, in agreement with what is expected in 
the AF basis based on a^. Note that to get the answer 
correctly one needs completeness. Had the second-order 
calculation been done in a truncated space of low energy 
orbitals (say near the Fermi energy) a nonexistent state- 
dependent shift would have been obtained to order (y) 2 - 
We mention this since we will perform some calculations 
with such truncations and will need to be careful. 

All the above was for the case when the spin-orbit 
coupling was constant over the sample. However, as 
pointed out by Brouwer, Cremers, and Halperin 6 , since 
the Rashba term depends on the confining electric field 
felt by the 2DEG, which in turn depends on the den- 
sity, one can fabricate additional gates over the sample 
to make the spin-orbit coupling inhomogeneous. In this 
case, the spin-orbit matrix element has a typical size 



*If one denotes a state by \ip) and its time-reversed state 
by \ip T ), and if one denotes the time-reversed operator corre- 
sponding to O by O t , then we have {tpi \O\1p2) = (1P2 \O t \ipl), 
and (V>i|6i6 2 |-!/>2) = {ip2\dZOZ\i>l). Also note that r T = r, 
and p T = —p. 
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(x)V#> an( l ^ e system will crossover directly into the 
fully symplectic ensemble at the energy scale ( ^) 2 g5. 

Let us now examine the effect of adding a — JS 2 term, 
which is shown to be the correct starting point in the 
appendix. The result is best seen in the AF basis. If 
we add a — JS 2 term to the hamiltonian of Eq. (4) in 
the original orthogonal basis, following Eq. (13) we end 
up with —J(Saf + SSaf) 2 where SSaf is the leading 
correction due to the action of U. The leading term 
— J(Saf) 2 commutes with a^ F but the corrections do 
not. They will typically produce effects of order (j) 2 ^, 
(say in relaxation of Sf F ) but the effects can be linear 
in (x)"^7f at an d near degeneracies in many-body states 
as will become clear soon. 

The effect of J may also be seen in the orthogonal 
basis, but not so easily. If we turn on a J and resort 
to many-body perturbation theory to order (y) 2 , J will 
enter energy denominators, in a calculation similar to 
the one leading to Eqn. (19). If we set J = in any 
of these results, we must recover the free-particle case. 
Expanding the denominators in powers of J/(e^ — e w ) 
will generate nontrivial effects proportional to powers of 
J. These will match the results described above in the 
AF basis. 

We now turn to our main purpose. 



B. Very low temperatures 

We will begin by defining the states of interest in the 
extremely low energy sector, T << T < a << S, where 
r is a typical single-particle level width due to coupling 
to the leads. This will allow us to get a firm grasp of 
all the relevant effects. Later we will incorporate the 
effects of T ~ a << S. Let us first turn off H so and 
Ez in Eqn. (1). Consider the A particle state. Figure 
1 shows the lowest energy state with (S, S x ) — (0, 0). It 
has two particles in an orbital <po at the Fermi energy 
(set to zero). The filled sea underneath can be ignored. 
The energy of this state is 



E o = 2 * + — A 2 - * J. 



(21) 



We could also form an S — 1 state, symmetric in spin 
and antisymmetric between the orbital 4>Q an( i t ne one 
just above it with a spacing 5, labelled <j>\. The energy 
of these states is 



E s =i, s=c =±i,o = S+^N 2 -2J. 



(22) 



Note that S — 2 J, the energy difference between singlet 
and triplet, can be much smaller than S. This is why 
we keep the orbital <j)\ but not higher ones. For now, 
when we do not have the Zeeman term, we will ignore 
the triplet state, assuming 5 — 2 J >> T. 



f 



(1/2,1/2) (1/2,-1/2) 
(0,0) (1,1) 

FIG. 1. Schematic of obitals and four of the low en- 
ergy states (labelled by (S,S X )) that play a key role at 
r<<T<<a<<5. The arrows denote spins. 



1 





Now for the A + 1 particle state. The extra particle 
will go to the orbital state <p\ with spin up or down, the 
two states being degenerate. The energy of these states 
is 



E 



5+| (A 



1) 



(23) 



By varying a gate voltage V g we can make the A and 
A + 1 particle states degenerate so the system conducts. 
Assume this has been done. Let us calculate the conduc- 
tance of this peak, using simple ideas. Later we turn to a 
master formula, based on rate equations 20,21 , that covers 
any situation where T >> T. 

We want an electron to jump in from the left with s x - 
independent width T^i and hop out of the right with 

U > 2 

s^-independent width r^i. The subscript on T labels 
the S values of the A and A + 1 particle states degen- 
erate at this peak. The dot has to be in the singlet N- 
particle state (and not in either of the doublet N+l- par- 
ticle states of the same energy) to receive this electron. 
Since all three states are degenerate at this peak, the 
probability for this is 1/3. The reduced conductance G 
(which has a trivial temperature dependence removed) is 
given by 



G 



2hkT 



-G = - 



2 L o,i 



A 



(24) 



(The logic behind this definition of G will be apparent 
later when we do a full-blown calculation of conduc- 
tance.) The subscript on G stands for the two states 
that are degenerate at this peak and the superscript re- 
minds us that the Zeeman term is absent: Ez — 0. The 
factor of 2 comes from sum over spins. The dependence 
on T's maybe understood as follows. The left and right 
contacts are like resistors in series while the T's represent 
their conductances. Thus we must add their reciprocals 
to obtain the reciprocal of the effective conductance. 
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Since wc want the electron to hop on the upper orbital 
described by <pi, 



T L f = T L/R \ML/R)\ U 



(25) 



where the arguments L/R refer to the coordinates of the 
two leads within the dot and T L / R represents an overall 
pre- factor for each contact. For now let us lump the 
overall scale T L / R into 4>i(L/R) 2 and write 



<-Tn 



= G° = 



2 \ML)\ 2 -\MR)\ 2 
3\ML)\ 2 + \MRW 



(26) 



Later we will return to Eqn. (25) that explicitly displays 
the tunnelling rate as a product of two factors, one that 
depends on the point contact and another that depends 
on the dot wave function. 

Let us now add a Zeeman term Ez and suitably 
adjust V g to follow this peak (which has been done 
experimentally 22 ). Once Ez >> T, only one of the spin 
states, with S r = +i is allowed in the dot. The reduced 



-\ is allowed in the dot 
conductance now becomes 



i \ML)\ 2 -\MR)f 



2\ml)\ 2 + \MRW 



(27) 



where the superscript on G signifies that the Zeeman 
term has a value strong enough to essentially filter just 
one of the spin-^ states, but not so large the spin-1 state 
(1,1) is competitive with (0,0). Compared to G°, the 
pre-factor | has become \ since we have just one spin 
channel (not two as earlier) and two degenerate many- 
body states (not three as earlier). 

The reduced conductance G has dropped to a value 
equal to 3/4 of the initial maximum: 



G 1 = -G°. 

4 



(28) 



The situation is sketched in Figure (2). 

Let us crank up the Zeeman term further till the spin- 
1 state (1,1) meets the (5,5) state and we again have 
three degenerate states. At this point denoted by a star, 
the conductance is given by 



G* 



\(G 1 + G 2 ) 



(29) 



where G2 is the conductance past this point where the 
spin-1 state (1, 1) is the ground state: 



G = Gi 



1 \ML)\ 2 ■ \MR)\ 2 
2\ml)\ 2 + \MR)\' 



(30) 



Note that since the extra fermion is placed in the lower 
orbital, only (f> enters G^- 

To get the above result we need only realize that at 
the triple degeneracy the two N-particle states each oc- 
cur with probability 1 /3 and the conductance via each is 
proportional to Gi or G2. 




FIG. 2. Schematic of G as function of Ez at zero spin-orbit 
coupling (q = 0) for very low T. Some rounding due to T > 
is assumed. The initial drop from G° is due to suppression 
of one spin orientation. The rise near E* z is due to the triple 
degeneracy. 



The behavior near E% depends on G 1 and G 2 . For 
instance if G 1 = G 2 = G, 



G * = 3 G - 



(31) 



In this case G first drops from G° to | of its value to G 1 , 
rises to G° and drops to G 2 . If however G 1 is less than 
G 2 /2, the rise will be monotonic with no bump. Likewise 
if G 2 is less than G 1 /2, the fall will be monotonic. Figure 
(2) describes the case G 1 = G 2 . 

For now let us focus on this case G 1 = G 2 , to get a 
feeling for our approach, which is to ask how the above 
mentioned behavior of G(Ez) is altered by a spin-orbit 
term. The general form of this term is 



7 A^y 
flu 



iB„ 



(32) 

where fi and v labels single-particle orbital states, \ 
and I the direction of spin along x and A^ lL , or B^ v are 
real antisymmetric matrices: 



A 1 = -A 



B T = -B. 



(33) 



The typical size of A^ u or £? M „ can be estimated to be 
f° r case °f h° m °g eneous spin-orbit coupling, 
and (^)Sy/g for the case of inhomogeneous spin-orbit 
coupling. For a typical small quantum dot of size 200nm 
with g 10, with the spin-orbit scattering length being 
2/ito, we find the spin-orbit matrix element to be typi- 
cally 0.03(5 for the case of homogeneous spin-orbit cou- 
pling, and 0.35 for inhomogeneous spin-orbit coupling 6 . 
These estimates should be considered crude, since the ex- 
pressions for various matrix elements have prefactors of 
order unity 5 which depend on the geometry. 
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Let us now focus on the very low energy sector spanned 
by (0,0), (1,1) and(i, i). Since H so conserves particle 
number, it does not mix the S = | states with S = 0, 1 
and neither does it mix the two S = | states with each 
other'''. In the S = 0,1 sector spanned by |0, 0) and |1, 1) 
the effective hamiltonian matrix assumes the form 



a 
a 



(34) 



where the single parameter a = A i (see Eq.(32)) char- 
acterizes the spin-orbit interaction. 

Note also that spin-rotation-invariant interactions be- 
yond the Universal Hamiltonian, such as Landau inter- 
actions, will appear in our basis as diagonal corrections 
to the already random energies. 

This operator clearly leads to an avoided crossing and 
selects a unique ground state in this sector. At the pu- 
tative degeneracy point the ground state with a^O, de- 
noted by | a, *) is unique and given by the antisymmetric 
combination: 



K*) 



[|0,0> - |1,1>] 



(35) 



The reduced conductance at very low temperature is 

i [\4>i (L)\ 2 + \ML)\ 2 ] [\MR)\ 2 + \MR)\ 2 ] 



G» 



4{\ML)\ 2 + \ML)\ 2 + \MR)\ 2 + \MR)\ 2 } 

(36) 



To derive this result we need to remember that there is 
now a probability of i for the dot to be in the N-particle 
state and that the widths are given, for example by 



= E!l<!< 11^(^)10,0) - <i, i|^(L)|i,i)| 2 

<7=± 1 



(37) 

and that in each of the two terms only one value of a 
contributes. 

If we are looking for the effects of a, we need to com- 
pare G with and without a. The answer depends on the 
wave functions via \4> /i(L/R)\ 2 . 

Again consider the simple case where all are wavefunc- 
tions are equal at both leads. It is readily seen that 



G*(a) = G 1 = G 2 , 



(38) 



i.e., there is no second bump at the would-be level cross- 
ing. This is an example of an order unity change in 



the conductance due to a small spin-orbit perturbation. 
However, one has to be at a temperature lower than a to 
see it. 

Since the first bump is there at Ez — 0, the absence of 
the second should be an unambiguous signal of a. Fur- 
thermore as T is raised beyond a (something we will con- 
sider in detail later) the bump should reappear (though 
rounded) since the symmetric combination of |0, 0) and 
|1, 1) (that got repelled upwards) also starts contributing 
to G. The requisite T will also tell us how large a is. 
While this is happening the bump at Ez = should be 
disappearing due to thermal effects. 

This happy line of thought ends when we realize that 
in general the four quantities \<Pq / i(L / R)\ 2 at the leads 
are not all equal and that the ratio G* (a) / G* is very sen- 
sitive to their values. We emphasize that even if we find 
experimentally that G 1 = G 2 we cannot say anything 
definite. For example G 1 — G 2 can also be realized if 



\MR)\ 2 = \ML)\ 2 = A 
\MR)\ 2 = \ML)\ 2 = B. 



(39) 
(40) 



In this case wc find 



Gi 



Gi 



G 1 



c;2 =-(VA/B + ^BjA) 2 >l. (41) 



Thus the ratio can be any number > 1 with large ra- 
tios resulting if either A or B is small. This behavior of 
Gi^ a can be understood by examining Eqn. (36). It is 
seen that if tunnelling amplitude is weak at one end and 
strong at the other for one wave function and the oppo- 
site is true for the other, Gi.„ is able to feed off the 

2 

bigger amplitude at both ends, unlike G 1 or G 2 which 
necessarily have a weak amplitude at one end. 

Thus if we are to say with confidence what effect a will 
have on the ratio G*(a)/G* at E*, we need to extract 
the four numbers \<p /i(L/R)\ 2 in a given dot on top of a 
given peak. Here is our proposal. 

Recall from Eqn. (25 )that 



T L { R = T L/R \ML/R)\' 

u, 2 

rff = t l/r \ml/R)\' 



(42) 
(43) 



where T l/ r are overall prefactors that control tunnelling 
rates at the left and right contacts and are determined by 
geometry, barrier heights and so on, and independent of 
the wave functions. Let us take the given dot and man- 
ually suppress tunnelling at the right contact by raising 
the barrier there, so that T R <<T L . The conductances 
assume the values G 1,2 (i?) given by 



fin the case S = § there would be an effect. Also, in metallic 
grains the situation is more complicated 10 



G\R) = -T R \MR)\ 2 
G 2 (R) = lv R \MR)\ 2 - 



(44) 
(45) 
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Their ratio gives us 



r = 



G\R) \MR)\ 2 



&{R) \UR)\ 2 ' 

Likewise we can set T L «T R and measure 



I 



G\L) \ML){< 



G*(L) \MLW 



(46) 



(47) 



Armed with these four numbers: (G 1 ,G 2 ,r,l) all 
measured at very low temperatures, we can solve for 
\^ n (L/R)\ 2 : 



l^iWI 2 

l<M^)l 2 

\UL)\ 2 

\MR)\ 2 



2G 1 G 2 {l-r) 

G 1 - rG 2 
2G 1 G 2 (r-Q 

G 1 - iG 2 
\ML)\ 2 
i 

\MR)\ 2 



(48) 
(49) 
(50) 
(51) 



While the right hand sides of the first two equations are 
not positive definite (and can be negative if G 1 and G 2 
are assigned arbitrary values) it can be shown that if we 
input any values for G that come from a real measure- 
ment, i.e., given by Eqs. (27,30) this will not happen. 
We will illustrate this point later. Note also our assump- 
tion that the wave functions in question are unaffected 
by what we do to alter the tunnelling rates. 

We can now express the conductances of interest in 
terms of knowns: 



G* = -i^ + G 2 ) 



G*(a) = 



G 1 G 2 (l+r)(l + l) 
2(G 1 + rlG 2 ) 



(52) 
(53) 



It is now easy to specialize to any case we want. If for 
example G 1 — G 2 = G and I = r = 1, 



G* = -G*(a). 



(54) 



Before concluding this analysis of the T « a case, we 
must ask what the (|, | ) level is doing in the meantime. 
In particular we do not want it to cross the (5,5) state 
before (0,0) crosses (1, 1), or better still before G settles 
down to G 2 . To avert this we need at least 



S' >S + J 



(55) 



where 8' is the energy gap between the orbital <j)\ and the 
one just above it, Since for dots with r s ~ 1 JJd ~ 
.3 15 , we would be fine in a system where say 5 ~ .75 and 
S' ~ 1 .3(5. If however the gap at the Fermi energy is larger 
than the one above it, we should expect the (5,5) — > 



(|, |) transition to precede the (0, 0) — > (1, 1) transition. 
The corresponding calculations are simply variants of the 
above. 

To test the picture described above we need to know 
where to look for putative degeneracies. Consider the 
peak associated with the (0, 0) — (5,5) degeneracy. The 
gate voltage V g (with leverage suppressed) at which this 
occurs is 



V, 



(0,0)-(±,±) 



(Ez) = Ei 1 — Eoq 

2 ' 2 

= l L(2N+l)- 



As a function of Ez it has a slope of 
hand when (1,1) wins we have 



Ez 
~2~' 



(56) 



2 • 



On the other 



(Ez) = E11 — Ei t \ 



= J L(2N+l)-5+\j+^-. (57) 

which has a slope of +5 as a function of Ez- 

At the point E* z = S — 2 J, the slope changes from — | 

to i, as shown in Figure (3) 19 . 

If we turn on a, we need to find eo(a), the ground 

state of the energy matrix in the (0, 0) — (1, 1) sector: 



a 

a S-2J-E Z 



(58) 



in terms of which the gate voltage to stay on the peak is 



K 



-<»> = f<2iV + l)-f 



Ez 



- e (a). (59) 



Due to the avoided crossing, eo (a) is lower by an amount 
a and V g is higher by an amount a at E* z . A sketch 
is given in Figure (3). Note that the results are most 
reliable near E* z where things depend on a. As we move 
away, results depend on a 2 and errors due to truncation 
of the orbital space will enter. These are expected to 
be small for 5 — 2 J <C 6, but we will not reproduce the 
arguments here. 



C. Higher Temperatures: r<<T~(5-2J<<5 

Now we consider a more detailed study of tempera- 
tures comparable to a but still considerably smaller than 
S. Thus we want to consider just six states: two at 
S = ^ and four at S = 0,1. In this larger space an- 
other parameter (3 = -B01, associated with the matrix 
B^ u of Eq. (32) enters. This matrix element mixes the 
states (S,S X ) = (0,0) and (S,S X ) = (1,0). Since (1,0) 
is split by about 6 — 2 J near the triple crossing, we can 
estimate its effects to be of order 1 /(S — 2 J) which is 
very small in the case of interest to us. We will therefore 
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V 



g 




FIG. 3. Schematic of V g as function of E z at T = 0. The 
straight and curved graphs correspond to a — and a ^ 
and approach each other asymptotically. The slope change at 
E* z represents the ground state transition (0, 0) — > (1, 1). 

initially ignore it in the analysis. Later we show a graph 
representative of its insignificance. 

As mentioned earlier, the cost of populating <pi with 
S = 1 is not 5 but 5 — 2J, and we assume T is such that 
many-body states involving 4>q and (pi are accessible but 
not those involving higher orbitals. 

The proper formalism for a finite T computation in- 
volves rate equations, developed by Beenakker 20 for a 
model with charging energy alone, and extended to the 
Universal Hamiltonian case by Alhassid, Rupp, Kamin- 
ski, and Glazman 21 . We follow their notation here. The 
formula for conductance is* 



G = 



hkT 



(JV+1) 



(AT) 



(60) 



where (i) i and j label many-body states with N and 
N + 1 particles with energies £j and q , (ii) = 

e< j N+1 ' — e j ) " e -F~eV s (where ep is the Fermi energy, cho- 
sen equal to zero and V g is the effective gate voltage) , (iii) 

P^ is the equilibrium probability of being in state i , 
(iv) fij = f(cij) is the Fermi function, with /(0) = \, (v) 
t]l is the fraction of the bias voltage dropped at the left 
contact with rjL + rja = 1, (vi) Tfj is the width for many- 
body transition i — ► j at the left contact when an electron 
jumps into the dot and (vii) 'F are functions defined by 



non-equilibrium probabilities in P^ and P, 



■>{N+1) 



to be 



in states i and j respectively in the presence of a bias 
voltage V 



P.- 



(N) 



(JV) 



I + eV(3^ 



(N) 



p(N+l) = p(N+l) 



1 + eVpV 



(N+l) 



(61) 
(62) 



"'"We differ with them on an overall sign for this formula. 



The values of 'F's are determined by demanding that the 
P's constitute a time-independent solution to the rate 
equations and their normalization condition. The result- 
ing equations are 



ij 



+ („ L r§ - m?*)] = o w (63) 



E^-^oRr 



ij 



•r£)(*f +13 



+ (rnXfj - mTy)] =0 Vj (64) 



£3 



(JV+1)^(JV+1) 



= 0. 



(65) 



Several comments are needed here 20,21 . While there 
seem to be i + j + 1 equations for the 'F's, one of the 
first i + j equations is redundant, which is where the last 
equation, from normalization, comes to the rescue. The 
result, Eqn.(60), which seems to depend on 77, in fact 
does not. Likewise the answer which seems to depend on 
the left lead in an asymmetric manner, is not. (There is 
an equally asymmetric expression in terms of the right 
contact which gives the same result 20 ' 21 .) All this is ex- 
emplified by the case where Eqns. (63, 64) hold for each 
term in the sum, i.e., 



vF 



(N) 



(JV+l) 



13 



Vi, j 



(66) 



in which case the expression Eqn. (60), for G reduces to 

,2 .... p L r^ 



G = 



__ ST pW f 



(67) 



Note that in the very-low-temperature case we consid- 
ered with just a total of three states, Eqns. (66 ) and 
(67) are applicable. § Our heuristic conductance calcu- 
lations (and the definition of G in terms of G) followed 
from above and the fact that /(0) = |. However, in our 
numerical calculations at temperatures comparable to a 
they do not hold and we solved Eqs. (63,64,65) numeri- 
cally. 

We now display the result of solving the above equa- 
tions for elevated temperatures using six states: two at 



§ In all cases we considered, either j had just one value and i 
had two or vice versa. Suppose j had just one value. For the 
one redundant equation we are allowed to ignore, we choose 
the middle equation (that sums over i). The first equation 
then gives Eqn. (66) for each i. 
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FIG. 4. Plot of G versus E z for the case G 1 = G 2 and 
I — 1.1, r = .9. Note that when a = (dotted lines), there 
is a bump at as T — > which broadens with T, while 
with a ^ 0, (solid lines, non-zero spin-orbit coupling) there is 
hardly a bump as T — > 0, but the bump grows with increasing 
T. 

5 = | and four in the S = 0, 1 sector. Since the graphs 
are functions of four variables, G 1 , G 2 , r, Z we need to de- 
cide how to sample them here. We have found that there 
are four main classes and give one from each.** 

In Figures (4, 5) we consider the case G 1 — G 2 . Re- 
call that in this case, in the absence of spin-orbit coupling 
(a = 0) there is a bump at E* z in the reduced conduc- 
tance G*. This bump will broaden with increasing T. 
When we turn on a there are two cases. If the ratios of 
Eqs. (46,47) satisfy r ~ I there will be no bump in G*(a) 
for T « a and as we cross T > a, the a/0 will look 
more and more like a = 0, as shown in Figure (4). If 
however r « 1 and I » 1, the case a/0 will have a 
bump at E* z even for T « a that will grow in height 
as T increases as shown in Figure (5). This bump will 
be robust even as the one at Ez = is smoothed out by 
temperature. 

Figures (6) and (7) consider the case of considerably 



FIG. 5. Plot of G versus E z for the case G 1 = G 2 and 
I = 4, r = \. Note that when a — (dotted lines) , there is 
a bump at E* z as T — > which broadens with T, while with 
a ^ 0, (solid lines) there is a prominent bump as T — > 0, 
which grows with increasing T. 

unequal G's. We choose for definitcness G 1 = 2G 2 . The 
first case with Z = 2.5 and r = 1.5 does not have a bump 
in G*, or G*(a) at low and high T, as shown in Figure 
(6). However if / are r are different enough from each 
other (we have chosen I = 4.0, r = 0.25) there is a robust 
bump in G* (a) that grows with increasing T as shown in 
Figure (7). 

In summary, the most dramatic manifestations of a 
appear when the r and I are quite different in size. It is 
also interesting that the effect of a is felt at values of T 
considerably higher than a. This is because, as explained 
in the paragraph below Eq. (41), both the states at the 
avoided level crossing can and do take advantage of large 
level widths. 

Finally we consider the effect of adding the term (3 in 
the spin-orbit interaction. A typical graph is found in 
Figure (8) where even a [3 twice as big as a is seen to 
make little difference. 



D. Experimental checklist 



**In choosing parameters for these graphs one cannot choose 
all the G 1 ,G 2 ,r and I arbitrarily. This is because even though 
every choice of the four |0| 2 will lead to positive values for 
G 1 ,€f 2 ,r and I, the converse is not true. Not every choice of 
G 1 ,Cf 2 ,r and I has an inverse image with all [</)| 2 ' positive . 
To keep the left hand sides of Eqns. (48,49) positive we must 
either satisfy I > r, and r < ~% < I, as in our choices above 
or a similar set of restriction for r < I. This is not a concern 
in extracting \<f>\ 2 from the experimental G 1 ,G 2 ,r and I, they 
will automatically generate the sensible positive underlying 
values of |0| 2 . 



Our paper began by looking for spin-orbit effects by 
starting with some input quantities and calculating ex- 
perimentally measured numbers in terms of them. How- 
ever, these input numbers, unlike say the electron mass 
or charge, involved wave functions at the leads, energy 
gaps and so on, and are not known a priori and the re- 
sponse of the system to turning on a is very sensitive to 
these quantities. 

We suggest a possible road map below to deal with 
this, being fully aware that our experimental colleagues 
will ultimately devise a better solution. 
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- a=0.00, T=0.015 

— a=0.02, T=0.015 

— a=0.00, T=0.05 

— a=0.02, T=0.05 




1.4 



1.2 



G i 



0.8 



0.6 





1 


i i | 




























- a=0.00, T=0.015 






— a=0.02, T=0.015 






-■■ a=0.00, T=0.05 






— a=0.02, T=0.05 


i i~ — ; — — 



0.1 0.2 0.3 0.4 0.5 



FIG. 6. Plot of G versus E z for the case G 1 = 2G 2 and FIG. 7. Plot of G versus E z for the case G 1 = 2G 2 and 
r = 1.5,1 — 2.5. There is no prominent bump with or without I = 4, r = \. There is no bump at a = (dotted lines), but 
a. there is one with a =fc (solid lines) which grows with T. 



• Go to the lowest temperature possible. 

• By studying the Ez dependence of V g 22 make sure 
the conductance peak describes the degeneracy of 
(0,0) and (|, |) which changes over to the degen- 
eracy of (1,1) and (i,|). Since a ^ 0, the rounded 
graph of Figure (3) will be seen. Since T « a, 
the rounding is not due to thermal effects.' 1 "'' 

• By retracing the asymptotically linear segments in 
V g locate E* z . 

• Since E* z = 8 — 2 J, find 5 if J is known or use its 
average value of .3(5. 

• Estimate the size of a by comparing the difference 
of V g 's between the two graphs at the cusp in Figure 
(3). 

• Measure G(Ez) and see if it looks like Figure (2). 
If so measure G 1 and G 2 . 

• By suppressing the overall tunnelling rate at the R 
(right) and L(left) ends, (with a controllable bar- 
rier) measure r and I, corresponding to the ratios 
of conductances. Bring tunnelling rates back to old 
value in the dot. 



^ A caveat here is that the parallel magnetic field also has 
the effect of enhancing the confinement of the 2DEG in the 
z direction 23 ' 22 , leading to a quadratic dependence in all gate 
voltages. This will have to subtracted before our suggestion 
can be implemented. 



• At E* z compare the predictions for G* (a) the value 
with a ^ to thejvalue_without spin-orbit coupling, 
namely G* = ^(G 1 +G 2 ). 

• Raise T and compare to theoretical predictions 
from the rate equations. All the parameters needed 
- \<p /i(R/ L)\ 2 , a, and 5 - are known by this point. 

• Try to raise the value of a by turning on a inho- 
mogeneities in the spin-orbit interaction using ad- 
ditional gates over the dot, following Brouwer et 
a/ 6 . 



III. DISCUSSION AND SUMMARY 

Spin-orbit couplings in quantum dots have been stud- 
ied extensively in the noninteracting limit in recent 
years 3-9 . In systems with homogeneous spin-orbit cou- 
plings, there is a new symmetry 4 which emerges at the 
Aleiner-Fal'ko (AF) scale e AF ~ (j) 4 gS, corresponding 
to a s z symmetry in a unitarily transformed basis 4 . This 
symmetry is broken at a scale {j) 6 gS. For small quan- 
tum dots with (j) ~ 0.1 and g ~ 10, one might therefore 
expect spin-orbit effects to be utterly negligible. 

Here we considered the problem when a robust ex- 
change term — JS 2 of the Universal Hamiltonian is 
present. There are order unity effects on the conductance 
at and near special degeneracy points, where two differ- 
ent many-body ground states with equal particle number 
but different total spin have the same energy. These ef- 
fects have unique signatures that can be produced only 
by the spin- orbit interaction. We have investigated the 
circumstances in which spin-orbit effects have the best 
prospect of making their presence felt, and proposed ways 
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1.1 




FIG. 8. Plot of G versus E z for the case G 1 = G 2 and 
I = 4, r = |. Both are at at a = 0.02 but one has (3 = and 
the other /3 = 2a — .04. Clearly /3 makes little difference. 

to measure the effect and compare it to theory. The ideal 
venue for studying the phenomenon discussed here is in 
dots where J has already brought the S = 1 state close 
to the 5 = state and one tunes through the degeneracy 
by a Zeeman coupling 22 . Our results are robust under 
the addition of non-universal Landau type interactions. 

Our calculations, done in the orthogonal basis, showed 
that in the presence of J spin-orbit terms enter the hamil- 
tonian to order (4)- We also saw how this phenomenon 
would appear in the AF basis, in which the spin-orbit in- 
teraction has been subsumed in the very choice of basis 
and a conserved quantity a AF emerges. In this basis the 
culprit is the — JS 2 term which does not commute with 
Haf, the non-interacting AF hamiltonian. Upon unitary 
transformation, one finds — JS 2 = - J(S AF + 8S AF ) 2 , 
where SS AF (which does not commute with a AF ) is or- 
der (t)/\/5 f° r homogeneous spin-orbit coupling. Wc 
have verified, by doing a calculation in the AF basis, 
that with the nonzero J, the AF degeneracy is indeed 
lifted and the spectrum coincides with what was found 
in the orthogonal basis. 

Although the spin-orbit term appears to a lower order 
in j in the presence of J than in the AF basis without 
J, its parametric dependence on the Thouless number 
g is very different for homogeneous spin-orbit couplings. 
Whereas in the AF basis a typical matrix element of the 
spin-orbit coupling is (j) 2 y/g, in the original basis its 
typical size is {k-)/y/§- The reason is evident from Eq. 
(14) 3 . In the original basis, the typical matrix elements 
of spin-orbit coupling averaged over the entire Thouless 
band are order ( j)^/<7, but the matrix elements between 
levels separated by 8 are of size I n tne AF 

basis, the spin-orbit matrix elements are those of (r x 
Paf) which does not suffer from this suppression. This 
is why Brouwer et af suggested that one could boost 



the spin-orbit term by making its coefficient 1/A space 
dependent. In this case the matrix elements of the spin- 
orbit coupling, even between neighboring states, would 
go parametrically as (j{)y/g~. 

Throughout this paper we have restricted ourselves 
to low temperatures T <C S, since that allows us to 
deal with only a few many-body states. Let us briefly 
consider the effect of raising T to be of order 5 or 
higher, and considering the effect on the conductance 
peak height distribution 13 . Excited states are populated 
at higher T, and these will suffer stronger spin-orbit ef- 
fects, since the typical matrix elements of the spin-orbit 
coupling increase with energy (Eq. (15)). In paral- 
lel with this increase, states of higher spin suffer larger 
spin-orbit effects, even holding the matrix element con- 
stant. It will be interesting to ask how this affects the 
peak height distribution at moderate temperatures, and 
in particular if it brings the theory into better accord 
with experiments 24 . 

We have concentrated here on small quantum dots in 
which e AF <C 5, but where {k)5/ \fg is still accessible at 
low temperatures. Our philosophy has been to keep g 
constant, and perturb in (4). 

For conceptual completeness L let us consider the other 
limit, when e AF = (j) 4 g<5 » 5. This regime would be 
accessible for very large dots. For noninteracting systems 
Aleiner, Fal'ko and co-workers 4 ' 5 have made predictions 
for conductance averages and fluctuations in open dots, 
which have been confirmed experimentally 25 . 

For closed dots with interacting electrons, one ex- 
pects that at energies very low compared to the AF 
scale e AF , the Universal Hamiltonian has only Ising ex- 
change interactions 11 , since only S AF is conserved (see 
appendix). As one increases the energy, one expects to 
see universal crossover effects in, for example, the spec- 
tral density of transverse spin excitations 12 . At energies 
much larger than e AF orthogonal symmetry is restored. 
Both these calculations started with the interaction term 
— J(S AF ) 2 , which, as we know, differs from the correct 
starting point (see appendix), — JS 2 . However, the dif- 
ference, which is at least first order in (5S AF , vanishes 
as g — > oo. In this limit, the previous analyses 11 ' 12 are 
correct. 

Finally, let us briefly touch upon the Stoner transi- 
tion in the presence of a small spin-orbit coupling. In 
the orthogonal class, one sees a series of metamagnetic 
transitions 14,15 with states of successively higher con- 
served spin as J — > 5, with the bulk Stoner transition 
being their accumulation point. In the presence of spin- 
orbit interactions, for any finite g these steps in S are 
rounded, and the transition will become a regular second- 
order quantum phase transition, where the order param- 
eter is not conserved (as we have seen above, even S AF 
is not conserved for finite g). On the other hand, if one 
works in the universal limit g — > co while keeping e AF /S 
fixed, there will be a sequence of metamagnetic transi- 
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tions between states with increasing conserved S z . For 
very large but finite g one can expect the steps of Sf F 
to be rounded and merge into a smooth curve as J — > 8. 
We hope to present a full analysis of the different regimes 
in a future publication. 

To conclude, we have found that in the presence of a 
robust exchange interaction J, one can, by tuning the 
Zeeman coupling 22 , readily expose order- unity effects of 
the spin-orbit coupling on the conductance at and near 
a many-body degeneracy. The effects have signatures 
that allow one to trace them back to the spin-orbit cou- 
pling. To access these effects experimentally one needs 
to reduce the temperature to be of the order of the spin- 
orbit coupling strength. Our results are unaffected by a 
modest amount of spin-conserving non-universal interac- 
tions. We have offered a proposal for a sequence of exper- 
iments to determine all the relevant parameters, that is, 
the nearest neighbor level spacing and the single-particle 
wave functions at both leads, which should enable exper- 
imentalists to test the predictions of this paper. In typ- 
ical small dots with 8 ~ O.lmeV, (-y) ~ 0.1 and g ~ 10, 
we hope that measurements to test the theory will be 
achievable in the near future. 
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original disordered basis with both the two-body Fermi 
liquid interaction and in addition a spin-orbit coupling. 
Without the spin-orbit coupling, this has been done by 
one of the authors and Harsh Mathur 16 . The results are: 
(i) Any two-body interaction which is a conserved quan- 
tity squared will not flow. This includes the charging 
interaction and the — JS 2 Stoner interaction in the Uni- 
versal Hamiltonian. (ii) Any other two-body interaction 
will flow, and flow to zero at weak coupling. This means 
that one can neglect the non-s-wave Landau interactions 
since they are irrelevant, and allows us to start with the 
universal hamiltonian — JS 2 coupling at the Thouless en- 
ergy. This argument also shows that higher channel spin- 
rotation-invariant Landau interactions do not affect the 
results of this paper, as long as they are weak. 

Coming now to the added spin-orbit coupling, this will 
grow in RG as the cutoff scale decreases, indicating that 
it is relevant for low-energy physics. At some cutoff scale 
the size of the spin-orbit coupling will rival that of the 
cutoff. This scale should be identified as caf- Below this 
scale, one cannot treat the spin-orbit coupling perturba- 
tively, but should start with exact states which incorpo- 
rate it. In the present paper, we are working in the case 
caf <C 8, allowing us to treat the spin-orbit coupling 
perturbatively. 

For the case caf 3> 8, neglecting the terms of order 
(^) 2 €af which break this S Zi af symmetry, there is a 
separate universal hamiltonian regime at energies much 
smaller than caf- This regime was identified by Alhassid 
and Rupp 11 and has only that part of the interaction 
which commutes with the kinetic term ~J Z S Z AF . One 
can obtain the same result by RG, this time diagonalizing 
the kinetic term with spin-orbit, and realizing that since 
the in-plane parts of the original Stoner interaction arc 
not conserved, they do flow, and flow moreover to zero 
leaving behind the —J z S zAF . 



Here we show that adding the — JS 2 is the correct 
effective theory at the Thouless scale Et as long as 
€af <C Et, and discuss the possible low-energy theories 
for energies much less than caf- 

Our reasoning is based on renormalization group 
(RG), which is an unbiased method for revealing the 
relevant and irrelevant terms, and is as follows: As is 
well-known, fcrmionic RG 26 in the clean system leads 
to Landau Fermi liquid theory 27 on a scale -Eflt much 
smaller than the Fermi energy. Since Et can be made as 
small as we want by increasing the size L, we can always 
achieve the condition Et <C Eflt- Therefore the cor- 
rect starting point at the Thouless scale is a Fermi liquid 
theory. This theory is parameterized by Landau interac- 
tion parameters 27 u m m — 0, 1, • • • in the spin singlet and 
triplet channels. For m = the Landau parameters are 
(up to a factor of 8) none other than the charging and spin 
exchange energies. Now consider carrying out RG in the 
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